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ABSTRACT 

We present an analysis of a Chandra ACIS-I observation of the massive galaxy cluster A644. This 
cluster was previously classified as a cooling flow, but no radio emission has been detected from 
its cD galaxy. Outside the core (R ~ 75 kpc ~ 0.03r v i r ) the hot ICM has properties consistent 
with a (relaxed) cool-core cluster out to the largest radii investigated (R ~ 415 kpc ~ 0.14r v i r ). 
Over this region the gravitating mass profile is described well by a Navarro-Frenk- White profile with 
concentration parameter, c = 6.1 ± 1.2, and virial radius, r v ; r = 2.9 ± 0.4 Mpc. However, inside the 
core the temperature and entropy profiles reverse their inward radial decline and rise at the center; 
the inner temperature profile is inconsistent with a constant at the 2.3a level. Although the core 
region does not display X-ray cavities or filamentary structures characteristic of radio-loud, cool-core 
clusters, the peak of the X-ray emission is offset from that of the centroid of the global X-ray halo 
by w 60 kpc. The position of the cD galaxy lies approximately between the X-ray peak and centroid, 
further testifying to a merger origin for the properties of the X-ray emission in the core. We discuss 
the implications of A644 and the small number of radio-quiet, cool-core clusters for the AGN feedback 
paradigm to suppress cooling flows in clusters. 

Subject headings: X-rays: galaxies: clusters - galaxies: halos - galaxies: formation - cooling flows - 
galaxies: clusters: individual: A644 



1. INTRODUCTION 

XMM and Chandra grating and CCD observations of 
the cores of clusters previously believed to be harboring 
coolin g flows find very little gas co o ling: below T ~ 1 keV 
(e.g.. iPeterson et alJ 12001 T2003 iTamura et al] 120011 
[ David et all 1200 lUMolendi fc Pizzolatol 120011 IXu et all 
I2002t lEttori et alJl2002t iBuote et al.ll2003aj) . This ma- 
jor discovery has led to the widely discussed picture 
where feedback on the hot intracluster medium (ICM) 
from an AGN in the central cluster galaxy is respon- 
sible for stifling the cooling flow, though the details of 
the heating process remain elusive (e.g., for reviews see 
iMathews fc Brighent,l2003lFabiaiJ20 04). About 70% of 
clusters previously classified as cooling flo ws possess ra - 
dio emission from the central galaxy (e.g.. lBurnsl [l990'). 
Chandra observations indicate that such clusters also 
have X-ray cavities and filamentary structures in the X- 
ray images that are related to the radio emission (e.g., 
iFabian et al.l2000tlM~cNamara et al-KOOttlBlanton et all 
l200lUMazzotta et al.ll200l iBirzan et al.ll2004|h 

What about clusters indicated to be cooling flows 
from previous low-resolution X-ray observations that do 
not display radio properties typical of cooling flows? 
We have previously analyz ed Chandra observatio ns of 
two suc h systems: A2029 llLewis et alJl2002L 12003) and 
A2589 ijBuote fc Lewis! 120041) . A2589 has no observed 
radio emission associated with its cD galaxy; and, al- 
though the cD in A2029 does possess a Wide-Angle- 
Tail (WAT) radio sourc e, such sources are generally no t 
found in cooling flows ijOwen et al J 11984 Burns! Il990|) . 
though their morphology provides evidence for relative 
motions between the radio source and ICM fEilek et alJ 
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iBurns et alJ I2002D . Both of these clusters have 
high density cores with temperatures significantly lower 
than at larger radii. (Although the relatively low-quality 
Chandra data of A2589 are consistent with an isother- 
mal gas, a new XMM observation clearly indicates that 
the temperature profile increases with radius - Zappa- 
costa et al. 2005, in preparation.) These clusters also 
display no spectroscopic evidence for cooling flows and 
have unusually symmetric X-ray images without the pro- 
nounced surface brightness irregularities observed in the 
cores of "radio loud" cool-core clusters fe.g.. lBirzan et alJ 
2004). Finally, cool -core clusters also appear to preferen- 
tially display central metallicity enhancements, while flat 
metallicity profiles are characteris tic of non-cooling flow, 
generally disturbed, sy stems (e.g.. lDe Grandi et al.l2004l 
iBohringer et al.ll2004() . A2589 and A2029 both have pro- 
nounced central metallicity enhancements further attest- 
ing to their comparatively relaxed states. 

Observations suggest that cooling flows are destroyed 
by large turbulent motions in the ICM since major 
merging activity in clusters anti-correla tes with cool- 
ing mass-flow rate l)Buote fc Tsail 11996 1 ). a conclusion 
supported by theoretical studies (e.g.. |Roet^ig er et alJ 
119961 iNorman fc Brvanlll999t iFuTita et alJl2004D Like- 
wise, theoretical s tudies (e.g., lLoken et alJ 119951: 
iRoettiger et al.llT996l) suggest that WAT radio morphol- 
ogy may also be a signature of hot shocks in on-going 
cluster mergers. Since the X-ray images of A2029 and 
A2589 are very regular all the way into their cores 
(though see and yet there is no evidence for ICM 
cooling or radio emission typical of cool-core clusters, in 
the context of the AGN feedback paradigm we are view- 
ing these clusters at a special time after they have settled 
down from past mergers but very soon before the onset 
of AGN heating. 
Here we examine a Chandra observation of A644, an- 
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other "radio quiet" cluster that was previously classified 
as a cooling flow ba sed on analy sis of both X-ray imagin g 
and spectral data <)Edge et alJll992t iPeres et alj fl998). 
A644 is the only cooling flow cluster with large inferred 
M (> 1 00 Mm yr" 1 ). which was not detected in the radio 
study of lBurnsI 1(1990. P 6cm < 3 xlO 21 W Hz" 1 ). We wish 
to examine how its X-ray properties compare to the other 
cooling flow clusters with atypical radio properties and, 
in particular, whether it demonstrates a clearer connec- 
tion between past merger activity and the suppression 
of a cooling flow. The redshift of A644 (z = 0.0704) 
corresponds to an angular diameter distance of 277 Mpc 
and 1" = 1.34 kpc assuming f2 m = 0.3, Oa = 0.7, and 
H = 70 km s" 1 Mpc -1 . 

The paper is organized as follows. In [J21 w e present the 
observation and discuss the data reduction. The anal- 
ysis of the image and spectral data of the ACIS-I are 
presented in ^ and respectively. We calculate the 
gravitating matter distribution in Finally, in Sj^l wo 
present our conclusions. 

2. OBSERVATIONS AND DATA ANALYSIS 

A644 was observed with the ACIS-I CCD array for 
approximately 30 ks during AO-2 as part of the Chan- 
dra Guest Observer program. We reduced the data us- 
ing the CIAO v3.1 (with CALDB v2.28) and HEA- 
SOFT v5.3.1 software packages. We followed the stan- 
dard CIAO threads 3 and reprocessed the level-1 events 
data to make use of the latest calibration information, in- 
cluding correcting for charge-transfer inefficiency and a 
time-dependent gain shift. From regions of least source 
contamination of the CCDs we extracted a light-curve 
(5.0-10.0 keV) to identify periods of high background. 
The observation was remarkably quiescent, and after re- 
moving a small period of modestly increased background 
rate, the total exposure used for subsequent analysis was 
28.9 ks. 

Since the cluster emission fills the entire ACIS-I field, 
we adopted for our default analysis the standard back- 
ground events files provided in the CALDB collected 
from suitable blank fields. For comparison to results ob- 
tained using these standard background templates, we 
also modeled the background directly. Using regions as 
far away from the cluster center as possible, we extracted 
a spectrum (with point sources masked out) and fitted 
it with a combine d source and background model (e.g., 
iBuote et al.l 12004') . For the source, we adopted a thin 
thermal plasma component (apec) with T = 6 keV and 
half-solar metallicity. The background was modeled with 
components for the Cosmic X-ray Background (CXB) 
and instrumental background: two soft APEC thermal 
plasmas and a hard power law for the CXB; a bro- 
ken power law and two gaussians for the instrumental 
background. Comparison of results using the alternative 
background is presented in ^1 

3. IMAGING ANALYSIS 

In Figure ^ we display the ACIS-I image of A644 in 
the 0.3-7.0 keV band. The image was corrected for ex- 
posure variations using an exposure map created for a 
monochromatic energy 1 keV, approximately the counts- 
weighted average energy of the spectrum. For pre- 

3 http:/ /asc. harvard.edu/ciao/threads 



sentation purposes only, the image was also adaptively 
smoothed using the CIAO task CSMOOTH with default 
parameter settings. The diffuse emission is seen to fill 
the entire ACIS-I field. The morphology of most of the 
cluster is quite regular and approximately elliptical in 
shape. 

Visual inspection of the central R ~ 30" region of 
the (smoothed) X-ray surface brightness reveals some ir- 
regular, non-azimuthally symmetric features (Figure^. 
These deviations from a symmetric configuration, though 
statistically significant, do not represent large fluctua- 
tions (< 10%) in the mean surface brightness level. In 
particular, the data do not provide clear evidence for 
distinct subclumps in the core region. 

We also searched for azimuthal fluctuations in the core 
by constructing a hardness ratio map (and also through 
spectral fitting in sectors, see jQ). We did not find any 
corresponding azimuthal fluctuations either in the core or 
on larger scales. Instead, we found the spectrally softer 
core gives way at larger radius to a spectrally harder 
surrounding medium. This radial variation is consistent 
with the radially increasing temperature profile (see JIJ ■ 

We measure the centroid and the ellipticity of the X- 
ray surfac e brightness using t he m oment method de- 
scribed by iCarter fc Metcalfe! £l980) and implemented 
in o ur previous X-ray studies of galaxies and clusters 
fe.g. . IBuote fc Canizaresll994l) . This iterative method is 
equivalent to computing the (two-dimensional) principal 
moments of inertia within an elliptical region. The el- 
lipticity is defined by the square root of the ratio of the 
principal moments, and the position angle is defined by 
the orientation of the larger principal moment. Follow- 
ing our previous s tudy of the elliptic ity of the Chandra 
data of NGC 720 l|Buote et al.ll2002fl we removed point 
sources and replaced them with smoothly distributed dif- 
fuse emission using the CIAO task DMFILTH. 

The centroid of the X-ray surface brightness shifts sig- 
nificantly (e.g., AR = 56 ± 6 kpc between 15-500 kpc), 
mostly directed toward the South. At the smallest radii 
the centroid is located within the region of image irreg- 
ularities in the core mentioned above. Outside of the 
core region the X-ray isophotes are very regular and 
moderately elliptical. The ellipticity gently falls from 
0.30±0.02 at a = 100 kpc to 0.243±0.004 at a = 500 kpc, 
where a is the semi-major axis of the elliptical aperture 
within which the moments are calculated. Over this ra- 
dius range the position angle (measured N through E) of 
the elliptical apertures is quite steady, varying between a 
maximum of 11 ± 2 degrees (at a = 100 kpc) and a min- 
imum of 6 ± 0.5 degrees (at a = 300 kpc). Overall, the 
image properties suggest a relaxed cluster outside about 
R w 100 kpc but a disturbed cluster within R w 50 kpc. 

We examined the radial surface brightness profile 
(with point sources masked out) located about the X- 
ray peak and also about the centroid computed within 
a = 500 kpc. In both cases the surface brightness profile 
is fitted fairly well out to R — 600 kpc by a (3 model: 
r c = 79 ± 2 kpc, f3 = 0.50 ± 0.01, X 2 = 200.4, 129 dof 
(peak) and r c = 118±3 kpc, (3 = 0.59±0.01, x 2 = 198.6, 
128 dof (centroid). The largest fit residuals (~ 20%) lie 
within the innermost region (r < 10") and are mod- 
est (< 5%) elsewhere. We note that the parameters 
r c and (3 increase systematically when we include data 
from increasingly larger radii or exclude data from the 
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Fig. 1. — (Left) Smoothed, exposure-corrected ACIS-I image of A644 in the 0.3-7.0 keV band before point sources have been removed. 
Overlaid contours are logarithmically spaced in intensity. (Right) Close-up of the R ss 45" region centered about the emission peak. 
The small circles indicate the (J2000) positions of the peak (08hl7m25.793s,-07d30m31.69s) and centroid (08hl7m25.106s,-07d31ml3.14s) 
determined from the X-ray data and the optical position of the cD galaxy (08hl7m25.6s,-07d30m46s). Note that the peak is computed by 
centroiding on the 30 brightest pixels. 



very central region. Both of these issues likely account 
for the different values obtained from previous ROSAT 
studies w hich included dat a at la rger radii; e.g., using 
the PSPC lEttori fc Fabianl l|1999ft obtain r c = 157 kpc 
and (3 = 0.69. 

4. SPECTRAL ANALYSIS 

We focus our attention on the spectral properties of 
the hot gas obtained within a series of concentric circu- 
lar annuli centered either on the X-ray centroid or peak. 
Below we also briefly discuss attempts to constrain az- 
imuthal variations in the spectral properties. The widths 
of the circular annuli were chosen so that the tempera- 
tures were determined to similar precision in each radial 
bin. For the centroid case we list the annuli in Table 
Annuli 1-4 contain 10 4 background-subtracted counts 
while bins 5-9 contain twice that amount in the 1.5-7 keV 
band. We restrict our analysis to energies above 1.5 keV 
because of uncertainties arising from excess absorption 
above the Galactic value discussed below, and because 
the gas temperature lies well above 1.5 keV within all 
annuli. For each annulus counts-weighted response ma- 
trices were generated using the CIAO tasks mkwrmf 
and MKARF. 

We fitted the background-subtracted spectrum with an 
APEC thermal pla sma modified by Galacti c absorption 
(6.4 x 10 20 cm- 2 . lDickev & Lockmanlll99ff) to each an- 
nulus. The free parameters are temperature, iron abun- 
dance, and normalization (emission measure). (The fit- 
ted abundance is technically a metallicity since we fit 
the iron abundance as a free parameter and the abun- 
dances of all the other elements tied to iron in their 
solar ratios. However, iron dominates the metallicity.) 
The spectral fitting was performed with XSPEC (11.3. li, 
IArnaudll 996) using the \ 2 minimization method. Hence, 
we rebinned all spectral channels to have a minimum of 



20 counts per energy bin (necessary for the validity of 
the x 2 method) and a signal-to-noise ratio of at least 3. 
The sol ar abundances in XSPEC are taken to be those 
given bv lGrevesse fc Sauvall l)1998t) which use the correct 
"new" photospheric value for iron which agrees also with 
the value obtain ed from solar-system meteorites (e.g., 
IMc WilliamI 1 1 997|) . To estimate the statistical errors on 
the fitted parameters we simulated spectra for each an- 
nulus using the best-fitting models and fit the simulated 
spectra in exactly the same manner as done for the actual 
data. From 100 Monte Carlo simulations we compute 
the standard deviation for each free parameter which wc 
quote as the "lcr" error. (We note that these la error 
estimates generally agree very well with those obtained 
using the standard A\ 2 approach in xspec.) All quoted 
errors are la unless stated otherwise. 

The radial profiles of temperature and iron abundance 
are shown in Figure Table Ogives the parameters ob- 
tained from the spectral fits for the centroid case; Fig- 
ure 01 shows the ACIS-I spectra of two representative 
annuli and the associated best-fitting one-temperature 
model. The simple one-temperature model provides a 
good description of the spectral data in all annuli. For 
R > 75 kpc the spectral fits for both the centroid and 
peak case agree very well. Over this region the tem- 
perature is consistent with a constant value ~ 8 keV 
as fo und also by previous studies of A6 44 with ASCA 
(e.g.. lWhitell2000HBauer fc Sarazinll2000D . In contrast to 
these previous ASCA studies, the Chandra data clearly 
reveal that the iron abundance declines with radius from 
a value of ~ 0.8Z Q at R ~ 75 kpc to ~ 0.25Z© at 
R ~ 400 kpc. Previous ASCA studies of A644 were 
consistent with a constant iron abu ndance profile (e.g., 
lWhitejl2000t iBauer & Sarazinll2000D within the large er- 
rors. We note also that the average iron abundance from 
ASCA (« Q.5Zq, scaled to our solar reference) obtained 
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r (arcsec) r (arcsec) 

Fig. 2. — {Left) Radial temperature and (Right) metallicity profiles computed within circular annuli oriented about two different center 
positions: peak (circles: blue) and centroid (squares: red), as explained in the text. The top axis indicates radial units in kpc. The centroid 
results are also presented in Tabic HI 

TABLE 1 

Parameters from the Spectral Fits 









R 


out 


T 




norm 




Annulus 


(arcsec) (kpc) 


(arcsec 


) (kpc) 


(keV) 


(solar) 


(10 -3 cm -5 ) 


(xVdof) 


1 


0.0 


0.0 


29.3 


39.3 


7.8 ±0.7 


0.53 ± 0.16 


2.54 ±0.07 


127.8/140 


2 


29.3 


39.3 


43.5 


58.5 


7.5 ±0.7 


0.78 ± 0.16 


2.63 ± 0.06 


118.3/146 


3 


43.5 


58.5 


56.6 


76.0 


8.1 ±0.6 


0.79 ± 0.16 


2.86 ± 0.09 


159.2/159 


4 


56.6 


76.0 


69.1 


92.9 


8.1 ±0.7 


0.83 ± 0.17 


2.94 ±0.10 


168.6/165 


5 


69.1 


92.9 


95.9 


128.9 


8.2 ±0.5 


0.74 ±0.12 


6.23 ± 0.13 


262.7/249 


6 


95.9 


128.9 


127.2 


170.9 


7.6 ±0.4 


0.55 ± 0.11 


6.01 ± 0.12 


235.9/225 


7 


127.2 


170.9 


180.8 


242.9 


8.5 ±0.4 


0.45 ± 0.09 


7.92 ± 0.16 


245.9/265 


8 


180.8 


242.9 


248.2 


333.5 


8.4 ±0.3 


0.40 ± 0.09 


6.78 ± 0.15 


254.4/248 


9 


248.2 


333.5 


364.1 


489.2 


8.1 ±0.5 


0.23 ± 0.09 


6.65 ± 0.10 


217.6/229 



Note. — Results of fitting a single APEC plasma emission model modified by Galactic absorption directly to the annular spectra (i.e., 
without spectral deprojection) over 1.5-7 keV for the centroid case. The norm parameter is the emission measure of the APEC model as 
defined in XSPEC: 10 — 14 (/ n e n p dV) / 4ir D 2 (1 + z) 2 with units cm -5 . The quoted errors are lcr computed using the Monte Carlo procedure 
described in El 



bv lWhitel l|2000D agrees well with our Chandra results. 

For R > 75 kpc the declining radial profile of iron 
abundance, approximately isothermal gas profile, and 
regular image morphology are ver y characteristic of re- 
laxed "cool core" clusters (e.g., iMolendi fc Pizzolatol 
l200lUDe Grandi et aljf2004[) . Inspection of Figure [3 re- 
veals that the region R < 75 kpc is more complex. The 
iron abundance profile is approximately constant for the 
centroid case, with marginal evidence for a central dip. 
In the peak case, the profile declines and then appears 
to increase in the central bin: the value in the central 
bin (1.15 ± 0.31 solar) is inconsistent with that in the 
adjacent bin at the rs 2.5tr level. 

The temperature profile for R < 75 kpc does not 
correspond to a typical cool-core cluster. In the cen- 
troid case, the temperature is consistent with isother- 
mal. In the peak case, the temperature declines inward 
to R ~ 50 kpc like a cool-core cluster but then rises 
significantly in the central bin (T = 8.9 ± 1.2 keV). The 
99% confidence lower limit on the temperature estimated 
using A\ 2 = 6.63 is T — 6.7 keV, still larger than the 



minimum value of T = 6.5 keV reached near R — 50 kpc. 
A power law fit to the temperature profile of the inner 
three radial bins is inconsistent with a constant value 
at the 2.3cr level and is highly inconsistent with the de- 
clining temperature central profiles characteristic of cool 
core clusters. 

Despite the good quality of the fits, we examined 
whether they could be improved further. Allowing abun- 
dances other than iron to vary separately did not yield 
noticeable improvement. Adding another temperature 
component also had no effect on the fits, even when en- 
ergies down to 0.5 keV were included. Consequently, 
we found no need for an additional multiphase cooling 
flow component in the spectra. Previous ASCA studies 
were divided be tween indicating a mas sive cooling flow 
(> 60Mq vr -1 . l^uCTj^Sarazinl |2000|) and no cooling 
flow (< 4OM vr-\ IWhitdl27K ! )0r . 

Using a standard model of a gas cooling at constant 
pressure as implem ented in our previous studies (e.g., 
iBuote et al.ll200 3a). our Chandra results indicate a best- 
fitting cooling rate of zero with 90% upper limits of 
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w 5 M© yr 1 in all annuli except annulus 1 (referring 
to the centroid case), where M = 7(< 11) Mq yr" 1 
(90% confidence). (We note that sometimes evidence for 
a cooling flow and multiphase gas in clusters has been 
claimed based on the analysis of CCD spectra extracted 
from a large aperture c overing, e.g., the enti re putative 
cooling flow region fe.g- lClarke fc et alJ200|) . However, 
the distribution of temperatures produced by the radial 
variation of a single-phase medium within the putative 
cooling flow aperture, as well as that contributed by the 
projection of higher temperature gas outside the aper- 
ture if the data are not deprojected, generally cannot be 
distinguished f rom a two-temperature medium or a cool- 
ing flow (e.g.. iBuotd 119991: IBuote et al.ll2003a|) . Hence, 
we have explored the need for multiphase gas using the 
narrowest annuli allowed by the data quality.) 

Although the hardness ratio map displayed no sig- 
nificant azimuthal variations O, we searched for az- 
imuthal spectral variation by dividing up the annuli into 
four equal sectors We rebinned the radial bins by a fac- 
tor of two to achieve enough signal for interesting con- 
straints. When fitting a single APEC component to each 
sector we found no evidence for azimuthal spectral vari- 
ations, consistent with the findings from the hardness 
ratio map. The temperatures and iron abundances in 
each sector within an annulus were found to be consis- 
tent within their « 90% confidence limits. In partic- 
ular, we do not find the > 15 keV temperature spike 
in a region a few arcminutes to the East of the center 
reported by iBauer fc Sara^ml IpOOOj) from ASCA. (The 
orientations of our sectors were oriented approximatel y 
to match those of sectors 2-5 of IBauer fc Saraziril2!M : | 

Finally, we have studied the sensitivity of the results to 
various sources of systematic error which we summarize 
below: 

Deprojection: We performed a non-parametric depro- 
jection using the "onion-peelin g" method as i mplemented 
in our previous studies (e.s.. IBuote) l2000|: IBuote et al.l 
l2003aHLewis et al.ll2003j) . That is, one begins by deter- 
mining the spectral model in the bounding annulus and 
then works inward by subtracting off the spectral contri- 
butions from the outer annuli. This technique tends to 



introduce fluctuations between parameters in adjacent 
radial bins. We found that to obtain results of com- 
parable quality to the projected analysis it was neces- 
sary either to fix the temperature at a constant value 
near 8 keV or restrict the radial variation in the the iron 
abundance so that its logarithmic slope was < ±1. Af- 
ter imposing either of these restrictions we obtain results 
consistent with those reported above for the projected 
analysis - both in terms of fit quality and parameter val- 
ues. (We note that these results for the temperatures 
and iron abundances obtained using our implementation 
of the "onion-peeling" method are very consistent with 
those obtained using the deprojection scheme provided 
in XSPEC using the PROJCT routine and require the same 
restrictions on the temperatures and iron abundances 
noted above.) 

Galactic Column Density and Bandwidth: When ener- 
gies below ~ 1 keV are included in the fits then a model 
assuming the nominal foreground Galactic column den- 
sity (6.4 x 10 20 cm" 2 ) overestimates the flux at low ener- 
gies, especially near 0.5-0.6 keV. If Ah is allowed to vary 
the fits are improved substantially within each annulus; 
e.g., in annulus 4 A% 2 = 45.9 for 228 dof. The best- 
fitting column densities are w 14 X 10 20 cm" 2 - about 
twice the Galactic value. Since the excess column densi- 
ties apparently do not vary with radius it is likely that the 
excess is associated with an underestimate of the Galactic 
column. (Using the background models described below 
we verified that the inferred excess column density is not 
the result of over-subtracting the soft background using 
the standard templates of nominally blank fields.) Our 
result is similar to that reported by an XMM study of 
A478 (|Pointecouteau et al.ll200l) . 

Compared to the results presented above for energies 
above 1.5 keV, the temperature and iron abundance are 
mostly unaffected when fitting down to 0.5 keV with a 
column density about twice the Galactic value. The most 
significant change is a systematic decrease of w 0.7 keV 
in the temperature values in all annuli, but the shape 
of the radial temperature profile is preserved. Although 
we expect the excess column density is of Galactic origin 
we preferred to confine our analysis to energies above 
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1.5 keV to avoid sensitivity to the precise value. 

Background: As described in for comparison to 
results obtained using the standard blank background 
fields we also constructed a model background by fit- 
ting the total (source plus background) spectrum of re- 
gions of the ACIS-I as far away from the center of 
A644 as possible. We found that results obtained us- 
ing this modeled background agreed extremely well with 
those of the standard blank fields. As expected, by far 
the largest differences were observed in annulus 9; e.g., 
T = 8.6 ± 0.6 keV using the modeled background com- 
pared to T = 8.1 ± 0.5 keV for the standard blank fields. 

Plasma Code: We investigated the sensitivity of our 
results to the plasma code using the mekal model. The 
quality of the fits and the temperature values were found 
to be very consistent within the la errors. The iron abun- 
dance values also were consistent with those obtained us- 
ing the APEC code, though the values obtained with the 
mekal code were generally smaller by 10% -20%. These 
results are consistent with prev ious studies l)Buote et al.l 
I2003albl: iHumphrev et ani2004jh 

Calibration: We explored the sensitivity of our results 
to the version of the Chandra calibration using previous 
versions of CIAO (v3.0.2) and the CALDB (v2.26). We 
found the results on the fitted spectral parameters to be 
consistent with those quoted above within the estimated 
« lcr statistical errors. 

5. GRAVITATING MASS AND GAS FRACTION 

To determine the gravitating mass distribution from 
X-ray observations requires the hot gas to be in hydro- 
static equilibrium. Both N-Body simulations and grav- 
itational lensing studies confirm the reliability of X-ray 
mass measurements, particularly for clusters with regular 
X-ray image morphologies (e.g., for a review see iBuotel 
|2004[) . The assumption of hydrostatic equilibrium should 
be valid outside the central ~ 75 kpc region of A644. 
Outside the core region the derived radial variations in 
the density and temperature of the hot gas agree for both 
the cases where the cluster center is defined as the surface 
brightness peak or centroid. Consequently, we exclude 
the central three annuli (defined in Table ^) from our 
default analysis, though for comparison we summarize 
results obtained using all annuli. 

The approach we use to calculate the gravitating 
mass distribution follows closely our p revious studies 
(jLewis et all 120031: iBuote k Lewis] 12004ft . As is stan- 
dard we assume spherical symmetry to provide a spheri- 
cally averaged mass profile appropriate for comparison 
to other observations and to the spherically averaged 
mass profiles obtained from cosmological simulations. In 
the equation of hydrostatic equilibrium we evaluate the 
derivatives of the three-dimensional gas density (p g ) and 
temperature (T) using simple parameterized models as 
discussed below. 

For our default analysis we projected parameterized 
models of the three-dimensional quantities, p g and T, 
and fitted these projected models to the results obtained 
from our analysis of the data projected on the sky (Table 
In this manner we obtained good constraints on the 
three-dimensional radial profiles of p g and T. Although 
our analysis closely follows our previous studies, we note 
two minor improvements. First, when fitting models to 
the radial profiles of the gas density and temperature to 



the data binned in annuli on the sky (as given in Table 
nj we integrated the models over each radial bin (rather 
than only evaluating at a single point within the bin) to 
provide a consistent comparison. Second, when project- 
ing models of the gas density and temperature along the 
line of sight we also included the radial variation in the 
plasma emissivity (A(T, Z-p c )). 

In Figure 0] we show the radial profile of the emission- 
weighted projection of p g 2 ; i.e., the "norm" parameter 
(see caption to Tabled of the APEC model divided by 
the area of the annulus. A single (3 model provides a 
good description of the data with fit residuals < 5% (Ta- 
ble EJ). Because the inner three annuli were excluded, 
a marginally larger value of (3 was obtained compared 
to the surface brightness fits (Q. Adding a second (3- 
model component provides small improvement in the fits 
with smaller residuals (< 3%) than the 1/3 model (Table 
0); the 2(3 model is displayed in Figure 0] Although the 
parameters for the 2(3 model are not as well constrained 
as those for the 1(3 model we adopt the 2(3 model for 
our analysis of the mass distribution because it describes 
the data a little better. Note that if the inner three an- 
nuli are included then the parameters obtained are con- 
sistent with those in Table El but are somewhat better 
constrained. For the case where the annuli are oriented 
about the emission peak the second component of a 2(3 
model (fitted to all annuli) matches very well the 1(3- 
model results obtained for the centroided case listed in 
Table El demonstrating the similarity of the gas density 
profiles obtained outside the core region for the peak and 
centroid case. 

The radial profile of the emission-weighted projection 
of T is shown in Figure 0] (middle panel) along with the 
best-fitting power-law fit. The simple power law is an 
excellent fit (Table EJ) and demonstrates that the tem- 
perature profile is consistent with isothermal within the 
estimated lcr errors. Notice that despite the inner three 
annuli being excluded from the fits the model describes 
those data points well. 

We constructed the mass profile using the models just 
described for the gas densi ty and temperatur e. Follow- 
ing previous studies (e.g.. iLewis et alJ 120031 . for each 
annulus listed in Table El we assign a single (three- 
dimensional) radius value, r = [(i?out + ^f/ 2 )/^ 2 ^ 3 \ 
where i?i n and i? ut are respectively the inner and outer 
radii of the (two-dimensional) annulus. The radial mass 
profile is plotted in Figure 0] using the 2/3-model for p g 
and the power law model for T. (We emphasize that al- 
though the mass values for the inner three data points 
are shown, being derived from the gas density and tem- 
perature models, they are not included in the fits.) 

The NFW model is a good smooth fit when exclud- 
ing the inner three radial bins (Figure 0} . We obtain 
X 2 = 6.9 for 4 dof with fit residuals < 15%. Because the 
mass data points are correlated it is necessary to recal- 
ibrate the x 2 statistic to examine goodness of fit using 
our Monte Carlo simulations O; i.e., for each set of pa- 
rameters for the gas density and temperature profiles ob- 
tained for each simulation, a mass profile is constructed 
and then fitted with an NFW model. These simulations 
indicate the the \ 2 value quoted above is within one stan- 
dard deviation of the expected value. 

In Table |31 we give the parameters of the NFW fit and 
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Radius (kpc) 

Fig. 4. — Models fitted to radial profiles (centroid case) as discussed in In each case the inner three data points are excluded from the 
fits. (Left) Chandra radial profile of the projected gas density squared (J p g 2 A(T, Zp c )di/A ann ) obtained by dividing the norm parameter 
of the APEC model (Table Q by the area of the annulus; A an n is the plasma emissivity evaluated using T and Zp c of the annulus. The 
diamonds represent the binned data and the solid line the best-fitting double-beta model. The fit residuals are plotted in the lower panel. 
(Center) Temperature profile of the binned data (diamonds) and best-fitting power-law model. (Right) Mass profile (circles and error bars) 
and best-fitting NFW model. 

TABLE 2 

Deprojected Radial Profiles of Gas Density and Temperature 



Gas Density 



r c PgQ ^c,2 Pgo,2 

Model (x 2 /d°f) (kpc) (3 (10 -26 g cm" 3 ) (kpc) /3 2 (1CT 26 g em" 3 ) 



1/3 12.1/3 120 ±8 0.64 ±0.02 2.27 ±0.12 

2/3 4.3/0 213 ±33 2.0 ± 0.3 2.19 ±0.24 352 ± 115 1.1 ± 0.3 0.82 ±0.17 



Temperature 



Tioo 

Model (x 2 /dof) (kcV) p 

power law 2.7/4 7.9 ± 0.4 0.04 ± 0.05 



Note. — These models were (emission-weighted) projected along the line of sight and fitted to the data about the centroid excluding 
the inner three annuli. Tioo is the temperature evaluated at r = 100 kpc and p is the exponent. Note that the best-fitting value of /3 for 
the 2/3 model was not well constrained and pegged at the upper limit we allowed in the fits. 



the estimated statistical and systematic errors. The esti- 
mates for the systematic errors are mostly based on our 
discussion in 21 an d we have also included the differ- 
ences obtained when using a 1/3 model to parameterize 
the p g profile. (Note we do not provide an estimated error 
arising from the deprojection method because the statis- 
tical errors are considerably larger when using the onion- 
peeling method.) The derived concentration parameter 
(« 6) and virial radius (« 2.9 Mpc) imply a virial mass 
of w 1.5 x 10 15 M Q appropriate for a massive galaxy clus- 
ter. The va lue of c is consistent with that expected fro m 
CDM (e. g ..lThoma,s 5^l2?M IfSiitsinirii et alJlSnol. 

We mention that within the central region (i.e., the 
central three annuli), where the assumption of hydro- 
static equilibrium is suspect, the mass density profile is 
flatter than NFW; i.e., the cumulative mass profile at 
small radii lies below NFW as seen in Figure 0] How- 
ever, for the case where the annuli are centered about 



the emission peak we find the mass profile within the 
central ~ 70 kpc is highly uncertain but consistent with 
NFW. (To obtain a reasonable smooth fit to the centrally 
peaked temperature profile in this case, we modeled the 
profile with two power laws with exponential cut-offs at 
both ends.) The smaller mass values in the core for the 
centroid case are attributed to the flatter density and 
temperature profiles obtained with respect to the peak 
case; i.e., when the annuli are not positioned about the 
emission peak radial differences in spectral quantities are 
smeared out. 

6. CONCLUSIONS 

Outside the core (R ~ 75 kpc ~ 0.03r vir ) the hot ICM 
has properties consistent with a (relaxed) cool-core clus- 
ter out to the largest radii investigated (R ~ 415 kpc 
~ 0.14r V i r ). The X-ray surface brightness is very regular, 
moderately elliptical in shape, with no evidence of sub- 



TABLE 3 

Parameters of NFW Fit to the Gravitating Matter Profile 



Best AStatistical ABackground APlasma ABandwidth ACalibration Ap g 

c 6.1 ±1.2 -0.4 +1.2 -0.4 +0.1 +0.4 
r vir (Mpc) 2.9 ±0.4 +03 ^03 ^(12 -0.1 -0.1 



Note. — The virial radius (r v j r ) i s defined to be t he radius where the enclosed average mass density equals 103p c (z) appropriate for the 
ACDM concordance cosmology (e.g., [E ke et all l998t). The inner three annuli are excluded from the fits. The "Best" column indicates the 
best-fitting value and "AStatistical" the la statistical error. The other columns refer to systematic differences in the best-fitting values 
obtained for different choices (see end of J3} in the background, plasma code, bandwidth, and calibration. For the bandwidth case the data 
are fitted down to 0.5 keV using a fixed foreground column density of 14 X 10 2t) cm~ 2 (see end of The final column refers to differences 
obtained when using a 1/3 model to parameterize the gas density. 



structure; previous ROSAT studies generally indicated a 
relaxed cluster on such scales (e.g.. lBuote fe Tsail ll996'). 
The radial profile of the density of the ICM is fitted well 
with a single /3 model, with a double /3 model providing 
slight improvement. The temperature profi le is consis- 
tent with isothermal or gent ly rising (e.g., I Allen et al.l 
l20niHVikhlinin fc et al]l2004l) . The entropy factor (S = 

k B T/n 2 e /3 , e.g., S = 797±21 keV cm 2 at R = 0.14r vir ) is 
consistent with those ob served in cool-core clusters (e.g., 
iPiffaretti fc et al.ll2004^ . The iron abundance declines 
with increasing radius (e.s. ADe Grandi et al.l feCKM'). The 
gravitating mass profile is well fitted by an NFW pro- 
file with reasonable concentration parameter for a halo 
of its virial mass. Finally, the gas fraction (0.076 ± 0.002 
at R = 0.14r v ; r ) is quite typical of relaxed clusters ob- 
serve d by Chandra an d produced in CDM simulations 
(e.g,. IAllen et ^ l200l . 

Inside the core, however, the Chandra data reveal an 
irregular, disturbed system. The X-ray image displays 
asymmetrical isophotes about the emission peak which 
are offset from the global emission centroid by « 50 kpc. 
When oriented about the emission peak, the temperature 
profile of the ICM declines towards the center (down to 
R ~ 50 kpc), very much resembling a cool-core cluster. 
But at small radii the temperature profile turns around 
and increases toward the center; i.e., the inner temper- 
ature profile is inconsistent with a constant at the 2.3tr 
level. Both the iron abundances and entropy profiles mir- 
ror the shape of the temperature profile with an excess 
at the center. 

Since A644 does not possess significant radio emission 
l)Burnslll99fl) , and the Chandra image does not show any 
X-ray cavities, the suppression of a cooling flow cannot be 
attributed to the present day feedback from an AGN. Nor 
is it easily ascribed to heating by a past radio outburst 
from the black hole in the cD, because the synchrotron 
lifetime of 1 GHz-emitting cosmic rays is ~ 10 s yr for 
a magnetic field strength of a few /zGauss in the cluster 
ICM. Since this energy loss timescale is comparable to 
the sound-crossing time for the cluster core (and thus to 
the time needed to heat the core) , the core should remain 
hot for about the same time as the 1 GHz radio emission 
remains visible after heating is terminated. It is even 
more difficult to explain heating of the "hot spot" (i.e., 
the central "hot" radial bin oriented about the emission 
peak - see Figure 2) by an AGN in the recent past, be- 
cause an even shorter time since the cessation of cosmic 
ray injection is required due to its smaller physical size. 



Nevertheless, a search for a steep-spectrum radio 
source would be worthwhile. The absence of low- 
frequen cy radio emission (e.g.. VLA Low-frequency Sky 
burvey; iCohen et alEOQ^ would be an even more power- 
ful test against the AGN heating model in A644, because 
100 MHz-emitting cosmic rays have three times longer 
lifetimes in comparable magnetic fields. Moreover, al- 
though the higher temperature in the "hot spot" appears 
to be accompanied by a sharp increase in iron abundance, 
supernovae are a most unlikely heating source for even 
this small portion of the core because the required num- 
ber of supernovae would need to be accompanied by a re- 
cent, galaxy-sized (> 10 9 M Q ) starburst. Neither a blue 
galaxy nor a luminous e mission-line nebu la are visible in 
the cluster core of A644 j|Hu et alJll98^1 . 

However, the morphological irregularities in the core 
of the X-ray image are very likely the result of merg- 
ing. The cD galaxy (2MASX J08172559 -0730455) is lo- 
cated ps 20 kpc South of the X-ray emission peak nearly 
along the direction of the centroid shift noted in It 
is also « 40 kpc displaced from the centroid of the X- 
ray emission computed from larger radii. Evidently the 
cD galaxy and core ICM are "sloshing" in the potential 
well of the main cluster, and this energy has prevented 
the establishment of a cooling flow. The metallicity peak 
measured within the innermost radial bin probably indi- 
cates that previous merging did not destroy the initial 
peak, created when it was spatially coincident with the 
cD galaxy. This is because, as argued above, the peak is 
unlikely the result of very recent supernova enrichment. 
Also, the visual appearance of the A644 galaxies is quite 
diffuse, stretching ~ 2 Mpc linearly across the sky, sug- 
gesting a system far from relaxation. This morphology 
strengthens the argument for an on-going merger. 

The prop erties of A644 (this paper). A2029 
llLewis et all I2IM 1200^) . and A2589 <)Buote fc Lewisl 
2004) have important similarities. Their X-ray images 
indicate a peak which is offset from the centroid of 
the global X-ray halo, with A644 displaying the most 
extreme shift (« 60 kpc offset) compared to the others 
(A2589: w 10 kpc, A2029: « 4 kpc). Other than 
these offsets, the X-ray emission is very regular in 
each case with no indications of cavities as are seen in 
cool-core clusters with AGN. 4 Two of these clusters 

4 We remark that [plarkc & ct al. (2004) have emphasized the 
existence of low-level (R3 10%) surface brightness deviations from 
a smooth elliptical model in the central ~ 10 kpc of A2029 using 
the AO-1 Chandra image, and claim that the WAT morphology 
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have no detected AGN radio emission, while A2029 
has a WAT morphology AGN, very unusual (unique 
to our knowledge) for a cool-core cluster. Finally, the 
lack of nebular emission lines in the cD galaxies are 
chara cteristic of non-cool-core clusters l|Gardiel et alJ 
1998). 

Of these three clusters only A644 possesses a central 
temperature peak. Although A644 appears to be unique 
in that its temperature profile declines toward the center 
and then rises, there are a few reports of clusters in the 
literature with temperature profiles measured by either 
Chandra or XMM that peak at their centers. The clus- 
ters A2218 HGovoni et a,l Jpll IPratt fc et al.ll20f)5T) and 
A3921 (Bclsolc ct al. 2005) are violent, advanced mergers 
with with pronounced asymmetrical surface brightness 
and temperature distributions uncharacteristic of A644, 
A2029, and A2589. Another clu ster with a reported 
centr al temperature peak is A401 ijSakelliou fc Ponmanl 
2004), which is in the earl y stages of a major merger 
(with A399). Nevertheless. iSakelliou fc Ponm an ( 2004) 
argue that the properties of A401 are determined primar- 
ily from its own hierarchical formation before interaction 
with A399. In this context it is interesting that A401 
has a weak radio halo since we have suggested that ra- 
dio halos form only in massive clusters where a violent 
merg er has proceeded fully into the cluster core ijBuotd 
2001). If the other clusters had a similar history, their 
lack of radio halos indicate they arc presently at much 
more evolved states than A401. 



In the paradigm wherein the energy output from an 
AGN disrupts cooling flows, it is expected that just prior 
to the onset of AGN feedback a cluster will have the char- 
acteristics of a cool-core cluster but without significant 
radio emission or X-ray cavities. Indeed, the clusters 
A644, A2029, and A2589 are excellent candidates for this 
pre-feedback state, with A644 being at a relatively early 
state and A2029 at the very latest evolutionary state, 
where an AGN outburst has already occurred but little 
or no heating from AGN-generated cosmic rays has com- 
menced. The extremely advanced evolutionary states of 
A2029 and A2589 make them ideal for studies of their 
core mass profiles using the assumption of hydrostatic 
equilibrium. 
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anitcorrelates with these small fluctuations. Using a more recent, 
deeper (80 ks) Chandra exposure we confirm « 10% deviations in 
the surface brightness (Zappacosta et al. 2005, in prep), but we 
do not confirm a strong anticorrelation with the WAT morphology 
as seen in other cool-core clusters. Compared to other cool-core 
clusters with obvious central disturbances these m 10% deviations 
in the surface brightness, which translate to even smaller (si 5%) 
deviations in the gas density, indicate that the core of A2029 is in 
a much more relaxed state. For comparison, we find that Perseus, 



a cool-core cluster with an obvious distrurbance in the core, dis- 
plays up to 50% fluctuations in the surface brightness. It should 
be added that, unlike A2029, Perseus also has strong azimuthal 
fluctuations in the temperature and metallicity which, however, 
apparently consipire to produce approximate azimuthal balance in 
the gas pressur e and therefore approximate hydrostatic equilibrium 
ISanders et al. 2004). 
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